APPROVED 
BY 

DRAFTSMAN 


O.G. FIG. 


CLASS 


SUBCLASS 




9 dyS36C& 



PCT/US99/11973 



1/5 



Silence detection 



Artifact detection 



Pre-detectors : Inspect the 
whole signal for non EIG 
areas. 



Pre-processing 



Transformation 

(Biased 
Autocorrelation) 




Sureness grade - How 
"voiced" is the current 
window ? 



L 



Pitch estimation 
( of the current 
window) 



/ 



Divide the signal into 
overlapping windows: 
Window length: 4 seconds. 

Overlap: 3 seconds 
Pitch estimation of the 

individual windows. 



Postprocessor 




Post processor: 
Processes all the estimates and 
'sureness" grades into a final 
estimate 




Output: 8 
Final Pitch Estimation / 



Fig. 1 
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PCT/US99/11973 



input: Raw recording from a Langerhans lsle£\ 
[EIG is the electro-physiological signal from the) 

islet. 



Fig. 2 



Amp 



Silence detector: 
For each second: 

n\dx(Sig) - m\n(Sig) 
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Yes 




Segment classified as^ 
"Silence". > 



Artifact detector. 
Signal compared to an adaptive 
threshold based on a 300 second 
recording. 
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-No 




Segment classified^ 
as "Artifact". , 



segment not classifies 
as an "Artifact" or a 
"Silence". 



Note: The exact thresholds vary greatly due to the measurement 
methods, the measurement equipment, and the conditions of the 
experiment. The thresholds must be fitted manually for each 
experimental set-up. 
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input: x(n) : A 4 second long analysis window 
N ( Window length ) = 400 
FS ( Sampling frequency ) = 100 ( Hz ). 
t = Index of the first sample in window. 
Input values in window: x(l) to x( I + N - 1 ). 




1 






Preprocessing: 
Not used in the preferred embodiment. 








Windowing: A rectangular window is used. 
w(n) = I 0<n< N -\ 
0 Otherwise 
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Fig. 3 



Normalized Biased autocorrelation transform. 



#^=AM 



fc=0 



Fundamental frequency Decision algorithm 

( See detailed flow chart ): 
Output: - The Fundamental period size( In 



m 



Period samples ) 
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/ Output: 
No / Analysis window classified as "Unvoiced" 
~ \ No fundamental frequency is estimated. 
V Sureness = 0. 



Pt(rn) = 



Compute Sureness: m PeriQd is marked here as m. 
l 



a/-H 



Z [x(i+l)w(i)] 
i-0 



A/-H-I 

Z [;rO+/+mM./+m)] 

J=0 



A/-H-1 
2\ n=0 
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Outputs: 

Window is classified as "Voiced". 
Sureness and Fundamental period size values. 
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(Input: <f> t (m) A 
Autocorrelation coefficients, j 
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Find the highest peak in the 
allowed fundamental frequency 
range. 

Allowed range: 0.3 (Hz) to 5 (Hz) 
Allowed range converted to 
samples: 20 (samples) - 333 ( 

samples). 
m Hi g h.stp.ak = max(#(m) m=20 ) 
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Find the 5 highest extreme points 
in the allowed range: 
m1,m2, m3 f m4, m5. 



(ml = 
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Fig. 4 
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Yes 




Output: 



'HightstPeak 
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Yes- 




Output: 
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Inputs: 

For each analysis window: Fundamental Period size, 
Voiced / Unvoiced classification, Sureness Index. 
Global: Detected Non EIG areas: Silences and artifacts. 




Convert the all Fundamental Period values, into Fundamental 
Frequency ( "Pitch" ) values. 
Pitch = 1 / Fundamental Period length ( In seconds ). 



Loop once over ail the analysis windows: 
Assign: Pitch = 0 ( Hz ), In all windows containing only silences. 
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Loop once over all the analysis windows: 
Assign: Window is "Unvoiced", In all windows containing 
Pitch > 4 (Hz) 
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Loop once over all the analysis windows: 
Assign Pitch values to all "Unvoiced" windows ( Windows 
containing artifacts are also treated as "Unvoiced" ). 
Pitch values are derived from a smart average, of the pitch values 
of the neighbouring windows. 



Loop once over ail the analysis windows: 
If Pitch of the current window, is half or double the average of the 
neighbouring windows, than Pitch = smart average of the 

neighbours. 
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Loop once over all the analysis windows: 
If Pitch of the current window, is 35% bigger or smaller than the 
average of the neighbouring windows, than Pitch = smart average 

of the neighbours. 



Output: 

A Fundamental Frequency ( "Pitch" ) value for 
each analysis window. 
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Fig. 5 



